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1 Introduction 

The Lie derivative is an extremely familiar operation in differential geometry. Given a smooth 
vector field V and tensor field 4>, the Lie derivative Cv<j> is certainly an intrinsic or invariant 
construction. In practise, we can take this to mean that the result is independent of any local 
coordinate formula that may be used. Alternatively, we may employ an arbitrary torsion-free 
affine connection V to write explicit formulae such as 

(t>bc ^ V^afoc + (V b V a )cP ac + (V c V a )ct> ba , (1) 

checking that the result is independent of the connection used. In this formula, the indices are 
'abstract indices' in the sense of Penrose [18]. It follows that the expression is coordinate- free 
whilst its freedom from choice of connection follows immediately from the formulae 

v a v c = v a v c + r ab c v b and v a <p bc = v a av - v ab d <t> dc - r ac Ubd, (2) 

capturing all possible torsion-free connections via the choice of tensor T ab c = T ba c . For further 
details on this point of view, see [15] . 

Writing the Lie derivative as Cy(f> suggests that we are regarding it as a linear differential 
operator acting on the tensor field <j) for a fixed choice of vector field V. However, it is clear 
from expressions such as ([1]) that we may equally well fix (f> and regard the result as a linear 
differential operator acting on vector fields V. Sometimes, this is the natural viewpoint. For 
example, if 4> ab is symmetric then we may rewrite ([T|) as 

{£ V <t>)bc = <t^bcVaV C + <\>ac^bV C + (V 1 c c\> ab )V C 

and now if <j) ab is a metric and we choose V a to be its Levi-Civita connection, then 

£ v cb = V a V b + V b V a = 0, 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http://www.emis.de/journals/SIGMA/symmetry2007.html| 
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which is the usual way of viewing the Killing equation on a Riemannian manifold. Of course, 
the best way of viewing Cv<j> is as a bilinear differential operator or as a differential pairing. It 
is invariantly defined on any manifold and is an example of the subject matter of this article. 

One can easily imagine other expressions along the lines of the right hand side of (H|) that 
turn out to be independent of choice of connection. Restricting the choice of connection leads 
to a greater range of expressions that might be invariant. In the next two sections we shall find 
examples that are invariant when the connections are deemed to be Levi-Civita connections for 
metrics taken from a fixed conformal class. 



2 First order symmetries 

Suppose L : E — > F is a linear differential operator between vector bundles E and F. A linear 
differential operator T> : E — > E is said to be a symmetry of L if and only if 

LT> = <5L for some linear differential operator 5 : F — > F. 

When L is the Laplacian on M 3 (and the bundles E and F are trivial) , its first order symmetries 
are well-known [T7]. Specifically, they comprise an 11-dimensional vector space spanned by the 
following 

, df df df df 
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dz dy dy dz dx 

f„(z*-x z -i?)?£ + 2zx?£ + 2zy?£ + zf, (3) 
dx dz dz dx dy 

where x, y, z are the usual Euclidean coordinates on ]R 3 . Following |10j . it is convenient to 
rewrite the general first order symmetry as 

f^V a V a f + l(V a V a )f + Cf, 

where V a is an arbitrary vector field of the form 

V a = -s a - m a b x b + \x a + r b x b x a - \x b x h r a (4) 

and C is an arbitrary constant. In fact, vector fields of the form @ on W 1 are precisely the 
conformal Killing fields, i.e. the solutions of the equation 

V a V b + V b V a = %g ab V c V c , (5) 

where g ab is the usual Euclidean metric on M. n . More generally, on M n the first order symmetries 
of the Laplacian may be written as / i— > T>yf + Cf, where 

V v f = V a V a f + ^(V a V a )f (6) 

for an arbitrary conformal Killing field V a . This expression certainly resembles the Lie derivative 
and we now make this precise as follows. If (f>bc- de is an n-form, then 



= (n + lXV^)^...^] = (V a V a )<j) bc ... de - [(W b V a )<t> ac ... de + ■■■ + (V e V a )<j) bc ... da \ , 
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where we are following [18] in employing square brackets to denote the result of skewing over 
the indices they enclose and also noting that a skew tensor with n + 1 indices on M n necessarily 
vanishes. Hence, the Lie derivative on n-forms simplifies: 

CvAc-de = V a X/ ' a <f>bc-de + (V 'bV^ac-de + ■■■ + (V ' e V a )cp bc ... da 

= v a v a fo c ... de + {V a V a )4> bc -de- 

Thus, if h is a section of the bundle A" of n-forms, then 

C v h = V a V a h + (V a V a )h. 

Therefore, the expression © coincides with the Lie derivative if / is interpreted as a section of the 
line bundle (A™)'"" 2 ^ 2 " (supposing that the manifold is orientable so that such fractional powers 
are allowed). In particular, we deduce that [Dy,T>y/\ = T^[v,w]i as * s a ^ so readily verified by 
direct computation. Also we conclude that, with this interpretation, the differential pairing T>yf 
between V a and / is invariant in the sense discussed for Lie derivative in Section [TJ 



3 Conformal geometry 

An alternative viewpoint on these matters is as follows. If V a is the Levi-Civita connection for 
a metric g ab , then the Levi-Civita connection V a for the conformally related metric g ab = ^t 2 g a b 
is given on 1-forms by 

Va^ft = V a (j) b - T a 4>b ~ Tb4>a + ^ C 4>c9ab, (?) 

where T a = V a logf2 and T c = g^T^- Thus, associated with a conformal class of metrics we 
have a restricted supply of torsion-free connections related by formulae such as ([2]) where 

T ab c = T a 5 b c + T b 5 a c - g ab T c (8) 

for closed 1-forms T a . It is convenient introduce a line bundle L on a general conformal manifold 
as follows. A choice of metric g ab in the conformal class trivialises L, allowing us to view a 
section of L as a function /. We decree, however, that viewing the same section by means of 
the trivialisation due to g ab = ^l 2 g a b gives us the function / = 0/. More generally, sections of 
the bundle L w transform by / = $l w f and are called conformal densities of conformal weight w. 
Since 

Va? = v a (n w f) = n w (v a f + wT a f) 

we may view L w as equipped with a family of connections related by 

Va/ = V a f + U,T a f. (9) 

On an oriented conformal n-manifold, the line bundle A n is also trivialised by a choice of metric. 
Specifically, we may use the inverse metric g ab to raise indices and then normalise the volume 
form e bc ...de so that e bc ' de e bc ---de = n - Since g ab = £l~ 2 g ab it follows that e bc ...d e = Q n (-bc---de and 
hence that we may identify A n = L~ n . It is easily confirmed that ([7]) implies 



V a <t>bc---de = V a (j) bc ...de ~ nT a C 



'be- -de 



for n-forms 4>bc—de, which is consistent with ([9]) for w = —n, as it should be. Combining ([7|) 
and © gives us a formula for the change of connection on A 1 <g> L w , namely 

V a <t>b = V a 4>b + (w - l)T a (f>b ~ ~£b4>a + ^ C 4>c9ab- 
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Let us now consider how the Laplacian A = V a V a acting on densities of weight w is affected 
by a conformal rescaling of the metric. We compute: 

V a V fe / = V a V b f + (w-l)T a V b f-T b V a f + (T c V c f)g ab 

and so 

Af = r 6 V a V fe / = n~ 2 g ab (V a V b f + (w- l)T a V&/ - T 6 V a / + (T c V c f)g ab ) 
= n~ 2 (V a V a f + (n + w- 2)T a V a f) 

= fT 2 (V a (V a / + wT a f) + (n + w- 2)T a (V a / + wT a f)). 
Here, we are using g ab to raise indices on the right hand side. Regrouping as 

Af = Q~ 2 (Af + (n + 2w - 2)T a V a f + w(V a T a + (n + w- 2)T a T a )f) , 
we see that when w = 1 — n/2 there are no first order derivatives in / and 

Af = n~ 2 (Af - (n/2 - l)(V a T a + (n/2 - l)T a T a )/). (10) 
On the other hand, the Riemann curvature tensor transforms by 

R-abcd = ®? (Robed ~ ^aedbd + ^bedad ~ ^bddac + ^addbc) 

where 

-afo = V a T b - T a T fe + ^T c T c g ab . 
It follows immediately that the scalar curvature R = g ac g cd R a bcd transforms by 

R = n~ 2 (R - 2(n - l)(V a T a + (n/2 - l)T a T a )). (11) 
From (110p and (jlll) we conclude that 

A/ - j^Rf = (A/ - j^Rf) ■ 

If we also absorb the factor of £2~ 2 into a change of conformal weight, then we conclude that 
Y = A — ^"l 2 !) R 1S a conformally invariant differential operator Y : L l ~ n l 2 — > L _1_ri//2 . This is 
the conformal Laplacian or Yaraahe operator. The peculiar multiple of the scalar curvature that 
must be added to the Laplacian to achieve conformal invariance is often referred to a curvature 
correction term. 

Now we are in a position to compare with the first order symmetries of the Laplacian on M n 
found in Section Of course, the curvature correction terms are absent for the flat metric 
on M. n . Therefore, there is no difference between the Laplacian and the Yamabe operator 
provided we restrict our attention to conformal rescalings that take the flat metric to another 
flat metric. But flat-to-flat rescalings are precisely what are provided by the conformal Killing 
fields ((5|). Specifically, the conformal Killing fields are those vector fields whose flows preserve 
the conformal structure on M. n and, moreover, the flat-to-flat conformal factors generated in this 
way are general [2]. The conformal invariance of the Yamabe operator on a general conformal 
manifold therefore implies the invariance of the Laplacian on M. n under Lie derivative Cy for 
conformal Killing fields V. Here, as we have just seen, £y should be interpreted as acting on 
the bundles L 1_n / 2 and L~ l ~ n / 2 . At the end of Section[2]we found that the symmetry T>y given 
by © may be regarded as the Lie derivative on (A n )( n-2 )/ 2n . These viewpoints now coincide 
because 

ji^n J_ l ~ n \. (A n )^ n_ 2 )/ 2n ^1— n/2 

The same reasoning applies more generally. The geometric interpretation of conformal Killing 
fields combines with the conformal invariance of the Yamabe operator and we have proved: 
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Theorem 1. Suppose that V a is a conformal Killing field on a Riemannian manifold, i.e. 
a solution of the equation (JS|). Then T>y given by © is a symmetry of the Yamabe operator. 
More precisely, 

where 

V v = V a V a + ^(V a V a ) and 5 V = V a V a + ^(V a V a ). 

In fact, it follows from proof of Theorem 1] that these are the only symmetries of the 
Yamabe operator. 

We have already observed that the differential pairing Dyf given by ([6]) is invariant in the 
very strong sense of being intrinsically defined on any manifold. Conformal invariance is a 
weaker statement but one that is easily verified directly. From ([2]) and j5]), 

V a V b = V a V b + T a V b - T b V a + T c V c S a b V a V a = V a V a + nT a V a 

whilst for / a density of weight w we have ([9]) . It follows immediately that 

V a V a f-^ a V a )f (12) 

is a conformally invariant differential pairing and for w = 1 — n/2 this agrees with ©. 



4 Higher order symmetries 

There are some second order symmetries of the Laplacian on 1™ that should be regarded as 
trivial. If T> = hA for some smooth function h, then AT) = 5A simply by taking 5f = A(hf). 
All symmetries preserve harmonic functions but T> = hA does this by dint of annihilating them. 
Ignoring these trivial symmetries, the second order symmetries of the Laplacian in R 3 were 
found by Boyer, Kalnins, and Miller [3j. They add an extra 35 dimensions to the 11-dimensional 
space ([3]). In |10j all higher symmetries were found in all dimensions and, ignoring the trivial 
ones, there is an additional finite-dimensional space of second order ones. More precisely, there 
are second order symmetries 

V v f = V ab V a V b f + ^(V a V ab )V b f + ^li$ +2) (VaV b V ab )f (13) 

for any trace-free symmetric tensor V ab satisfying 

V a V bc + V b V ca + V c V ab = ^ ( gab V d V cd + g bc V d V ad + g ca N d V bd ) . (14) 

Such V ab are called conformal Killing tensors (of valence 2) and on W 1 form a vector space of 
dimension 

(n-l)(n+2)(n+3)(n+4) 
12 

As detailed in [10], this is the 'additional space' of second order symmetries referred to 
above. As with first order symmetries ([6]), the pairing between valence 2 symmetric trace- free 
tensors V ab and densities / of weight 1 — n/2 given by the right hand side of (|13p is conformally 
invariant under flat-to-flat conformal rescalings. Unlike the first order case, however, the pair- 
ing (|13[) is not invariant under general conformal rescalings of a general Riemannian manifold 
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but, like the Laplacian, becomes so with the addition of suitable curvature correction terms. 
Specifically, for V ab any trace-free symmetric contravariant tensor and / of weight w 

V« b V a V b f - ^(V Q OV 6 / + I ^^(V V 6 0/ + ^^R ab V ab f, (15) 

where R ab is the Ricci tensor, is a conformally invariant differential pairing on an arbitrary 
Riemannian manifold. When w = 1 — n/2, we obtain 

Vyf = V°»V a V h f + ^r 2 {V a V ab )V b f + 4(^^+2) (VaV 6 V° b )/ - ^^6^/ 

as a conformally invariant analogue of the flat symmetry operator (|13p . Whether this operator 
actually provides a symmetry of the Yamabe operator in general, however, is unclear. By analogy 
with Theorem \T\ one might hope that if V ab were a conformal Killing tensor (I14p , then it would 
follow that 

( A " W^) R ) Vv = S V<A- ^R) (16) 

for T>v as above and 

Svf = V ab V a V b f + ^(V a V ab )V b f + ^(VaVftV*)/ - ^R ab V ab f. 

This is currently unknown. The problem is that there does not seem to be any useful geometric 
interpretation of V ab being a conformal Killing tensor in parallel to V a being a conformal Killing 
field. In principle, the validity or otherwise of (|16p should boil down to a calculation once the 
differential consequences of (|14|) are determined. Furthermore, the differential consequences 
of (|14p can all be found by prolonging the system, as described in [5] for example. Nevertheless, 
it will surely be a difficult calculation. 



5 Further examples of invariant pairings 

The simple conformally invariant pairings we have found so far, namely (|12p and (|15p . may be 
extended by allowing V a and V ab to have a general conformal weight as follows. It is convenient 
to set 

<&ab = — ^ {Rab ~ iRgab) ■ 

n — 2 n ' 

Proposition 1. Suppose that V a has conformal weight v and f has conformal weight w. Then 

(v + n)V a V a f - w{V a V a )f (17) 

is conformally invariant. Suppose V ab is trace-free symmetric with conformal weight v and f 
has conformal weight w. Then 

( n + v + 2 )( n + v + l)V ab V a V b f -2(w-l)(n + v + l){V a V ab )V b f 

+ w(w - l)(V a V b V ab )f + w(n + v + w){n + v + 2)<5> ab V ab f (18) 

is conformally invariant. 

Proof. These are easily verified using the formulae for conformal rescaling developed in Sec- 
tion El ■ 
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Notice that there are several special cases occurring for particular values of v and w. The 
form of (|17|) suggests setting w = or v = —n and, when we do, we find 

(v + n)V a V a f is invariant \/V a / h- ► V n / is invariant 

and 

— w(V a V a )f is invariant V/ ==> V a t—> V a V a is invariant. 

These are familiar invariant linear differential operators. The first one is the exterior derivative 
d : A — > A 1 and the second may be identified as d : A n_1 — ► A n . More interesting are the 
consequences of setting w = 1 and v = —n — 1 in (|18|) . We obtain conformally invariant linear 
differential operators 

/•-» V a V 6 /-i(A/) 5a6 + $ ab / and V ab ^ V a V b V ab + <5> ab V a \ respectively. 

The theory of conformally invariant linear differential operators is well understood [T| Efl 1121 H4| 
and includes these two as simple examples. Other special values of the weights v and w show 
up as zeroes of the coefficients in ()18f) . We have already seen that for w = there is an 
invariant differential operator / ^ V a /. Similarly, for v = —n — 2 the operator V ab i— ► V a V ab 
is conformally invariant. 

It is also interesting to note that, conversely, we can build the invariant pairing (I18p from 
these various linear invariant operators. Specifically, if V ab has weight v and / has weight 
then, at least where / does not vanish, 

f-(n+v+l)/wyab ^ _ n _ 1 &nd f -(n+v+2) /wyab ^ _ n _ 2 . 

Therefore, 

V a V b (f-( n+V+1)/w V ab ) + <S> ab f-(n+v+l)/wyab and Vb ^-(n+v+2)/w v ab^ 

are invariant. If the second of these is multiplied by f 2 / w we obtain a vector field of weight — n 
and conclude that 

V a (f 2/w \7 b (f~ {n+v+2)/w V ab )) 

is also invariant. It follows, therefore, that the combination 

w(n + v + 2)(n + u + ^J/^+i-H^)/™ (V a V fe (/- (n+1,+1)/w y afe ) + $ afe /-(™+t ) +i)/» y a6^ ) 
- u;(n + v + l)(n + v + 1 + u; )/(™+^+«')/«'(v a (/ 2 /«'v 6 (/- (n+1 ' +2)/w 'y afc ))) 

is invariant. A short calculation reveals that this expression agrees with (I18p and, although we 
supposed w 7^ and also that / ^ in carrying out this derivation, the final conclusion is valid 
regardless. These tricks reveal some sort of relationship between invariant differential pairings 
and invariant linear differential operators, at least when one of the quantities to be paired is 
a conformal density. 

These tricks are unavailable more generally but there are, nevertheless, many conformally in- 
variant differential pairings. For example, if V ab is symmetric trace-free with conformal weight v 
and <p a has conformal weight w, then 

(n + v + 2)V ab V a <j )b - (to - 2){V a V ab ) ( t> b . 

is invariant. Operators such as this one and (|17f) are referred to as first order, not only because 
there are no higher derivatives involved than first order but also because there are no cross terms 
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of the form W ix V</> where IX denotes some algebraic pairing (we shall come back to this point 
soon). The theory of first order invariant differential pairings on a conformal manifold, or more 
generally on a manifold with AHS structure, is reasonably well understood thanks to recent work 
of Kroeske [16]. Following the approach of Fegan p3] in the conformal case and Cap, Slovak, and 
Soucek [7] more generally, Kroeske classifies the first order invariant pairings provided that the 
weights of the tensors involved are so that invariant linear differential operators are excluded. 
This corresponds well to the relationship that seemed to be emerging in our examples. For 
higher order pairings, however, Kroeske [16] finds additional unexpected pairings and is able to 
present a satisfactory theory only in the flat projective setting. 

6 Some precise formulations 

When E and F are smooth vector bundles on some smooth manifold, it is well-known that linear 
differential operators D : E — > F of order < k are in 1-1 correspondence with vector bundle 
homomorphisms J k E — > F where J k E is the /c th jet bundle of E. Indeed, it is usual to abuse 
notation and also write D : J k E — > F, sometimes adopting this viewpoint as the definition of 
a linear differential operator. There are canonical short exact sequences of vector bundles 

-f Q k A 1 ® E -» J k E -» J k ~ x E -» 0, (19) 
where Q k A l is the bundle of symmetric covariant tensors of valence k and the composition 

Q k A 1 ®E^ J k E F 

is called the symbol of D. For further details, see [19] for example. 

Differential pairings may be formulated similarly. Suppose another vector bundle V is given 
and we wish to define what is a differential pairing E x F ^ V . The jet exact sequences (|19h 
may be written as 

J k E = E + A 1 & E + Q 2 A 1 (g)E + Q 3 A 1 ®E + ••• + Q k A 1 E, 

meaning that J k E is filtered with successive quotients as shown. It follows that 

2 A X E F 
A 1 (g> E <8) F 
J k E <g> J k F = E&F + + A 1 <g> E (8) A 1 ® F + ••• 

E A 1 <g> F 

E ® Q) 2 A l F 

and, following [16], we define the bi-jet bundle J k (E,F) as the quotient of J k E J k F corre- 
sponding to the first k + 1 columns of this composition series. There are bi-jet exact sequences 
starting with 

A 1 E F 
0^ © -► J\E,F) ->£0F->O 
£0 A 1 0F 

and 

Q 2 A 1 ®E®F 

© 

0^ A 1 ®E®A 1 ®F -> J 2 (E,F) -> J l (E,F) -> 0. 

© 

S O 2 ^ 1 i 7 
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A differential pairing E x F — > V of order k may now defined as a homomorphism of vector 
bundles J k (E,F) — > V. Notice that it is the total order of the operator that is constrained to 
be less than of equal to k. For first order pairings, for example, we are excluding terms of the 
form Vtp ix (written in the presence of chosen connections). For second order pairings we 
exclude V 2 ip ixi V0, Vtp M V 2 (/>, and V 2 tp ixi V 2 0. The symbol of a pairing is defined as the 
composition 

k 

0O i A 1 ®£®O H A 1 0i? — ► J*(£,F)— >V. 
3=0 

Several other basic notions for linear differential operators immediately carry over to bilinear 
differential pairings. Suppose, for example, that we consider homogeneous vector bundles on 
a homogeneous space G/P. If E and F are induced from P-modules E and F, then J k (E,F) 
will also be induced from a P- module, say J fc (E,F). If V is induced from V, a G-invariant 
differential pairing will correspond to a P- module homomorphism J fc (E, F) — > V. In particular, 
when G = SO(n+l, 1) and G/P is the conformal sphere (as detailed in [9] for example) then these 
P-module homomorphisms correspond to flat-to-flat invariant differential pairings. This gives 
an approach to classification as adopted in [16]. Nevertheless, it remains a difficult problem. 

7 Other sources of invariant pairings 

Firstly, there are other conformally invariant differential operators whose symmetries give rise 
to invariant differential pairings. The square of the Laplacian is treated in [11] and the Dirac 
operator in [13]. There is also another source of pairings derived from the structure of the 
symmetry algebra as follows. Let us consider the symmetries of the Laplacian as in Sections [2] 
and HI It is clear that the composition of symmetries is again a symmetry. In particular, if we 
compose two first order symmetries ([6]) on M n then we find [10] 

V v V w f = V v@w f + \v [vm f - j^^iV, W)f + ±v a w a Af 
where T>v@w is given by (fT3|) and 

(V ® W) ab = \V a W b + \V b W a - ±g ab V c W c , 

[V, Wf = v b v b v a - w b v b v a , 

(V, W) = (n + 2)(V b V a )(V a W b ) - ^(V a y a )(V^ b ) - n±l V * VaVbW b 

_ n±2 W a Vg y b yb + A(F a VT). 

The last of these is constant when V a and W a are Killing fields on W 1 and in this case coincides 
with their inner product under the Killing form in so (n+ 1, 1). Remarkably, there is a curvature 
corrected version 

(n + 2)(V b V a )(V a W b ) - ^(y a V a )(V b W b ) - n±l V * Va V b W b 

_ n±2 ^^^^b + A( y aW a) _ ^^^yay^b + RV a W a 

that provides an invariant differential pairing on arbitrary vector fields. 

This approach to constructing invariant pairings is sometimes called quantisation: the symbol 
of an operator is specified, one attempts to build an invariant operator with this symbol, and 
then one composes these operators. Further examples are given by Duval and Ovsienko [8] in 
the conformal case and Fox [15] in the projective case. 

Finally, there is a general construction of invariant pairings for very particular weights called 
cup products by Calderbank and Diemer [6] . This construction applies in any parabolic geometry 
but further discussion is beyond the scope of this article. 
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